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ABSTRACT 

The recent discovery of high frequency oscillations in giant flares from SGR 1806-20 and SGR 
1900+14 may be the first direct detection of vibrations in a neutron star crust. If this interpretation 
is correct it offers a novel means of testing the neutron star equation of state, crustal breaking strain, 
and magnetic held configuration. Using timing data from RHESSI, we have confirmed the detection 
of a 92.5 Hz Quasi-Periodic Oscillation (QPO) in the tail of the SGR 1806-20 giant flare. We also find 
another, stronger, QPO at higher energies, at 626.5 Hz. Both QPOs are visible only at particular (but 
different) rotational phases, implying an association with a specific area of the neutron star surface 
or magnetosphere. At lower frequencies we confirm the detection of an 18 Hz QPO, at the same 
rotational phase as the 92.5 Hz QPO, and report the additional presence of a broad 26 Hz QPO. We 
are however unable to make a robust confirmation of the presence of a 30 Hz QPO, despite higher 
countrates. We discuss our results in the light of neutron star vibration models. 

Subject headings: stars: magnetic fields—pulsars: individual (SGR 1806-20)—stars: neutron—stars: 
rotation—stars: oscillations—X-rays: stars 


1. INTRODUCTION 

The Soft Gamma Repeaters (SGRs), objects that ex¬ 
hibit recurrent bouts of gamma-ray flare activity, are 
thought to be magnetars - ne utron stars with mag¬ 
netic fields greater than 1 0 14 G iDjmggi^^^HigmggSS, 
1^22 ! iThomoson fc Duncanlll995t iWood^^ThomDSoiT 
2004 ). On rare occasions SGRs exhibit giant flares, 
hugely energetic events with peak fluxes in the range 
10 44 — 10 46 ergs s _1 . The giant flares consist of a short 
spectrally hard initial peak, followed by a softer de¬ 
caying tail lasting a few hundred seconds. Pulsations 
with periods of several seconds are visible in the tail, 
and reveal the neutron star spin period. Their pres¬ 
ence is thought to be due to a fireball of ejected plasma, 
trapped near the stella r surface b y the strong magnetic 
field IThomoson fc Duncanill995h . 


Three giant flares have been detected since the advent 
of satellite-borne hi gh energ y detectors: in 1979, from 


SGB 0526-66 (|M?y,^ e 

al. 

[1979|); in 1998, from SGR 

1900+14 (Hurlev et ^1. 

1999 

); and in 2004 from SGR 

1806-20 (Terasawa et al 

2005; Palmer et alj[200a). The 


catastrophic magnetic instability that powers the giant 
flares is thought to be associated w ith large-scal e frac- 


turin 

s of the neutron star crust ((Flowers & Ruderman 

1977 

Thonmson & Duncan 19951 200j Schwartz et al. 

2005 

. This will almost certainly excite global seismic 


vibrations iDuncanlll 998fl : terrestrial seismologists regu¬ 
larly observe such global modes after large earthquakes 
such as the 2004 Sumatra-Andaman event iPark et al . 

EMI. 

The 2004 flare from SGR 1806-20 was the most en¬ 
ergetic ever recorded. Analysis of da ta from the Ross i 
X-ray Timing Explorer (RXTE) by llsrael et~all l|2005lf 
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revealed a transient 92.5 Hz Quasi-Periodic Oscillation 
(QPO) in the tail of the flare, associated with a particular 
rotational phase. The presence of 18 and 30 Hz features 
was also suggeste d, although with lower significance. 
llsrael et~ahl (2005) suggested that the 30 Hz and 92.5 
Hz QPOs might be toroidal s h ear vibrations of th e crust 

E£^makgT)^^h(^^_^83_^cDgmipt2£l al. ] mu 

iStrohmaved 119911 lDuncanlll998t iPiro 11200511 . These 
modes are particularly promising in terms of both ease of 
excitation and coupli ng to the magnetic field to modulate 
the X-ray lightcurve iBlaes et al. 11198 9th 

Motivated by this work. iStrohmaver fc Wattsl l2005fl 
re-analysed RXTE data for the 1998 flare from SGR 
1900+14. QPOs were found at 28, 54, 84 and 155 
Hz, again associated with a particular rotational phase. 
These QPOs can plausibly be identified with a series of 
toroidal modes. The association with a rotational phase 
away from the main peak (where emission from the sur¬ 
face is obscured by the trapped fireball) implies that we 
are seeing emission from a particular area of the stellar 
surface or bundle of magnetic held footpoints. 

The frequencies of crustal modes such as the toroidal 
modes depend on neutron star mass and r adius, crusta l 
rigidity, and magnetic held configuration l|Dun can; 1998). 
For this reason the high frequency QPOs offer a power¬ 
ful new diagnostic of neutron star properties provided 
that an accurate mo de identihcation can be made (see 
IStrohmaver fc Watts SI for an initial attempt to 
constrain SGR parameters). 

In this paper we present a timing analysis of the SGR 
1806-20 hyperhare using data from the Ramaty High En¬ 
ergy Solar Spectroscopic Imager (RHESSI), and discuss 
our results in the light of the crustal vibration model. 

2. OBSERVATIONS AND DATA ANALYSIS 

RHESSI is a solar-pointing satellite whose primary ob¬ 
jective is the st udy of sola r hares in the energy range 3 
keV - 17 MeV iLin et al. 112002 1. The detection system 
comprises nine cyrogenically-cooled Germanium detec- 
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tors, divided into front and rear segments. For solar flares 
all direct photons under 100 keV should be stopped in 
the front segments. Most direct photons above 100 keV 
are captured in the higher volume rear segments. Lower 
energy photons reach the rear detectors only indirectly, 
by scattering (there is a Beryllium scatterer embedded in 
the front segments, intended for use in polarization stud¬ 
ies). The rear segments also capture albedo flux from the 
Earth. 

The SGR 1806-20 hyperfl are was ~ 5° off-axis for 
RHESSI llBoggs et al. ifiooH) . The detector saturated 
during t he p eak, but recorded the decaying tail in its 
entirety llHurlev et al. I l|2005P : Figure El . Although the 
flare was not directly in the RHESSI field of view, most 
photons in the front segments would have been direct. 
Given RHESSI’s native time resolution of 1 binary p,s 
(2~ 20 s) these events are clearly suitable for high fre¬ 
quency timing analysis. The rear segment flux, by con¬ 
trast, comprises scattered photons from the front seg¬ 
ments, direct photons entering through the walls of the 
spacecraft, and albedo flux. The latter, which could be 
as much as 40-50% of the direct flu x in the energy range 
of interest dMcGonnell et al.1l2004D . has a severe impact 
on timing analysis. At the time of the flare RHESSI was 
passing the limbs of the Earth (as viewed from the SGR). 
Albedo flux is limb-brig hten ed, part icularly if the incom¬ 
ing flux is polarized (IWillis et al. mu. This means 
that a large fraction of the detected photons could have 
incurred additional delays of up to ss 0.02 s, smearing out 
signals above ~ 50 Hz. Note that although countrates in 
the rear segments exceed those recorded by RXTE, coun¬ 
trates in the front segments are slightly lower. It should 
also be noted that scattering from the spacecraft walls 
and the Earth will cause the photon energies recorded 
by RHESSI, particularly in the rear segments, to deviate 
from the true energies of the incident photons. Quanti¬ 
fying this effect precisely is extremely difficult. For this 
reason we use broad energy bands in our analysis, and 
urge some care in interpreting the recorded photon ener¬ 
gies. 

We started by extracting event lists from the RHESSI 
data, excluding only events occurring in a 2s period 
« 270s after the peak of the flare when an attenu¬ 
ator is removed (the associated spike introduces spu¬ 
rious variability, particularly in the front segments). 
Timing analysis was carried out using the i? 2 statisti c 
jBuccher^et ak JlQSflllStrohmaver fo Markwardt 11999(1 . 
Ilsrael et al.l (l2005D showed that the presence of the high 
frequency signals was dependent on the phase of the 
7.6s rotational pulse; the signals appeared most strongly 
at phases away from the main peak. Similar phase- 
depende nce was also observed in the SGR 1900+14 hy¬ 
perflare (iStrohmaver fc Watts 112005 ). As such we have 
conducted a phase-dependent analysis. 

We searched for phase-dependent QPOs by folding 
data of a given rotational phase from N p pulses, gen¬ 
erating power spectra that are averaged to a frequency 
resolution Aza The distribution of noise powers is a y 2 
distribution with 2 N degrees of freedom, where N = 
N p AvPAQ, P is the rotational period and A*I> is the 
phase window under consideration (0 < A4> < 1). We 
searched over a range of A4>, N p . and energy bands for 
any signals with significance > 3 ct. 

We started by searching for signals in the range 50- 



Fig. 1.— Light curves. Top: Front segments, 25-100 keV band. 
Center: Rear segments, 25-100 keV band. Bottom: Front seg¬ 
ments, 100-200 keV band. The plots show the main peak and 
decaying tail with the 7.6s double-peaked pulse profile. The spike 
in the front segments at 270s is due to the removal of an attenuator. 
Zero time corresponds to 21:30 UTC on December 27 2004. 


1000 Hz, using only data from the front segments. In 
this range the noise profile is Poissonian. We find only 
two signals that meet our search criterion. 

The first, for photons with recorded energies in the 
range 25-10 0 keV, is the QP O at 92.5 Hz previously re¬ 
ported bv Ilsrael et all ( 2005 + shown in Figure [21 This 
signal, which we detect only at a rotational phase away 
from the main peak, is st rongest » 1 50 -^_260 s after the 
initial flare. As noted bv Ilsrael et al.i (: 2005 b this occurs 
in conjunction with an increase in unpulsed emission. At 
Av =1 Hz the QPO is resolved; at Av =2 Hz it is not. 
We estimate the significance of the Av =2 Hz power us¬ 
ing a x 2 distribution with 68 degrees of freedom, which is 
the distribution expected based on the number of inde¬ 
pendent frequency bins and pulses averaged. The peak 
at 93 Hz has a single trial probability of 2 x 10 - '. Apply¬ 
ing a correction for the number of frequency bins, inde¬ 
pendent time periods, and rotational phases searched we 
arrive at a significance of « 1 x 10~ 3 . That th i s is low er 
than the significance reported bv Ilsrael et alJ (2005) is 
to be expected, given that the RHESSI front segment 
countrate is lower than that of RXTE. A search for the 
signal in the RHESSI rear segments indicates that the 
signal has indeed been smeared out due to albedo flux. 
Fitting the QPO with a Lorentzian profile we find a cen¬ 
troid frequency 92.7 ± 0.1 Hz, with coherence value Q 
of 40. The integrated RMS fra ctional ampl itude is 10 ± 
0.3 %, in good agreement with Ilsrael et al.l : 2 005: . 

The independent detection with RHESSI of the 92.5 
Hz QPO is a strong confirmation of the RXTE findings. 
Using the significance quoted bv Ilsrael et alJ (' 2005 !. we 
can compute the probability of getting two apparent de¬ 
tections at the same frequency, time and phase, due to 
noise alone, given the number of trials. If we do this we 
find that the detection of the 92.5 Hz QPO has a com¬ 
bined significance of > 6cr, an extremely robust result. 

The second detection, for photons with recorded ener¬ 
gies in the range 100-200 keV, is of a narrower QPO at 
626.5 Hz (Figure^. This signal is strongest in the period 

















































3 




Fig. 2.— Average power spectra from 2.27 s intervals (0.3 cycles) 
centered on different rotational phases, computed using photons 
from the front segments with recorded energies in the range 25- 
100 keV. The upper curve was computed using 15 successive 2.27 s 
intervals, « 150 — 260 s after the main flare, at a rotational phase 
that includes the secondary peak and part of the DC phase. The 
frequency resolution is 1 Hz. The middle curve shows the same 
spectrum with 2 Hz frequency resolution. The QPO at 92.5 Hz is 
clearly visible. The lower curve is for the same time period but 
is an average of rotational phases zb 2.27 s away from the 92.5 Hz 
signal phase: no QPOs are detected. Characteristic error bars are 
shown for each spectrum. 


50-200 s after the main flare. As before, the signal is only 
detected at certain rotational phases; but in this case the 
rotational phase where the signal is strongest is centered 
on the leading edge of the main peak (Figure 0J. We esti¬ 
mate the significance of the peak using a y 2 distribution 
with 114 degrees of freedom, which is the distribution 
expected based on the number of independent frequency 
bins and pulses averaged. The peak at 626.5 Hz has a 
single trial probability of 7.7 x 10~ 9 . Applying a correc¬ 
tion for the number of frequency bins, independent time 
periods, and rotational phases searched we arrive at a 
probability of chance occurrence of « 6.6 x 10 -5 . Fitting 
the QPO with a Lorentzian profile we find a centroid fre¬ 
quency 626.46 ± 0.02 Hz, with coherence value Q of 790. 
The integrated RMS fractional amplitude is high, at 20 
± 3 %. 

Extending our search to lower frequencies, in the range 
1-50 Hz, we can now use data from both the front and 
rear segments. In this range the noise spectrum is not 
Poissonian. If we focus only on the time perio d 200-300s 
after the main flare, where fTsrael et, all ( 2005 1 reported 
possible QPOs at 18 and 30 Hz, we do not find features 
that meet our search criterion. Extending our search to 
earlier times, however, we do find evidence for low fre¬ 
quency features at the rotational phase where the 92.5 Hz 
QPO was found (Figure 0. At 1 Hz resolution we find 
broad features at ~ 18 and 26 Hz, with a slightly weaker 
feature at « 30 Hz. Assigning a significance to these pu¬ 
tative detections requires us to fit the noise powers. We 
fit a constant plus power law model, excluding the candi¬ 
date peaks to avoid biasing the fit. We then divide by the 
continuum model and multiply by 2 to obtain the nor¬ 
malised spectrum shown in the lower panel of Figure 0 
We compute significances using a 2 Hz resolution power 


Fig. 3.— Average power spectra from 2.27 s intervals (0.3 cycles) 
centered on different rotational phases, computed using photons 
from the front segments with recorded energies in the range 100- 
200 keV. The upper curve was computed using 19 successive 2.27 
s intervals, « 50 — 200 s after the main flare, for a rotational phase 
centered on the main peak. The frequency resolution is 1 Hz, and 
the QPO is clearly visible at 626.5 Hz. The lower curve is for the 
same time period but is an average of rotational phases zb 2.27 s 
away from the main peak: no QPOs are detected. Characteristic 
error bars are shown for each spectrum. 
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Fig. 4. — Folded pulse profiles (background subtracted), showing 
the rotational phase where the high frequency QPOs are detected. 
The top panel shows the average pulse profile for photons from the 
front segments, in the 25-100 keV band, for the period where the 
92.5 Hz QPO is detected most strongly (2 full cycles are shown). 
The lower panel shows the average pulse profile for photons from 
the front segments, in the 100-200 keV band, for the period when 
the 626.5 Hz QPO is detected most strongly. The hatched regions 
indicate the rotational phases where the QPOs are detected. 


spectrum (not shown), to take into account the breadth 
of the candidate QPOs. If we neglect the uncertainty 
associated with the normalisation, and assume that the 
noise is Poisson, the single trial probabilities of the pow¬ 
ers measured at 18 Hz and 26 Hz are 3.9 x 10~ 8 and 
2.7 x 10 -9 respectively. Accounting for the number of tri¬ 
als, the probabilities of chance occurrence are 8.8 x 10” 6 
for the 18 Hz QPO and 6.1 x 10”' for the 26 Hz QPO. The 
uncertainty associated with the noise model will reduce 
the quoted significances by some factor, but unless the 
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Fig. 5.— Average power spectra from 2.27 s intervals (0.3 cy¬ 
cles) centered on different rotational phases, computed using pho¬ 
tons from the front and rear segments with recorded energies in 
the range 25-100 keV. The upper curve was computed using 22 
successive 2.27 s intervals, « 60 — 230 s after the main flare, for 
the rotational phase where the 92.5 Hz QPO is observed. Several 
possible QPOs are apparent at low frequencies, and we show the 
best-fit power law plus constant model of the continuum. The cen¬ 
ter curve is for the same time period, but is an average of rotational 
phases zb 2.27 s away from the main peak: no QPOs are detected. 
In the lower panel we have renormalised the power spectrum using 
the continuum model. Broad QPOs at 18 Hz and 26 Hz are readily 
apparent, as is a weaker feature at 30 Hz. Note that there is no 
longer a distinct QPO peak at 92.5 Hz, as there was in Figure 2, 
although there is a broader feature at this frequency. The fact that 
the peak at 92.5 Hz has been smeared out illustrates the effect of 
albedo flux in the rear segments on high frequency signals. The 
frequency resolution in all panels is 1 Hz. Characteristic error bars 
are shown. 


noise mo del is seriously in er ror, t hese detections are ro- 
bust (see Ivan der Klis ( I1989T) and llsrael fc Stella! (Il996f) 
for further d iscussion of th is issue). An 18 Hz feature was 
seen bv [Israel et al„ ( 2005 ): it would be interesting to see 
if a closer analysis of the RXTE dataset reveals the 26 
Hz feature. Fitting a Lorentzian profile we find for the 
lower frequency QPO a centroid frequency 17.9 ± 0.1 Hz, 
Q of 10, and integrated RMS fractional amplitude 4.0 ± 
0.3 %. For the higher frequency QPO we find a centroid 
frequency 25.7 ± 0.1 Hz, Q of 9 and an integrated RMS 
fractional amplitude of 5.0 ± 0.3 %. 

The weaker 30 Hz feature does not meet the 3er crite¬ 
rion at either 1 Hz or 2 Hz resolution. As such we are 
unable to make a robust c onfirmation of the 30 Hz QPO 
suggested by llsrael et ~ahl (.2005), despite higher coun- 
trates. One might boost the significance of this feature 
by combining RXTE and RHESSI d at asets ; but since the 
time window reported bv llsrael et al. (1200511 is different, 
some caution is warranted. 

3. DISCUSSION 

The strong rotational-phase dependence of the QPOs 
in the tail of the SGR 1806-20 hyperflare provides com¬ 
pelling evidence that all of the oscillations are assocated 
with particular regions of the stellar surface or magneto¬ 
sphere. In this section we discuss the candidate rnecha- 
n isms i n more deta il. 

llsrael et all (l2005fl identified the 92.5 Hz QPO with 
the l = 7, n = 0 toroidal mode of the crust. This in- 


tepretation was supported by the apparent detection of a 
QPO at 30 Hz, the expected frequency of the fundamen- 
taH = 2 mode. Although we cannot confirm the detec¬ 
tion of the 30 Hz QPO, the interpretation of the 92.5 Hz 
QPO is probably secure. If the m agnetic fie ld of ~ 10 15 
G in ferred from timing analysis (IWoods Kr. Thompson! 
12004 1 is accurate, then for the 92.5 Hz QPO to be an 
l = 6 toroidal mode the stellar mass would have to be 
< 1 Mq for even the softest equation of state. To be 
an l = 8 mode, the neutron star mass would have to 
exceed that of SGR 1900+ 14 b y > 0.3M© (see the dis¬ 
cussion in iStrohmaver fc Watts I ( 2005 11. An alternative 
to the toroidal modes, with a frequency in the correct 
range, is the crustal interface mode, but its behavior in 
the pr esence^ o f str ong magnetic ^fields^ re quires further 

The detection of a QPO at 626.5 Hz is exciting, as 
there are several crustal modes that could have frequen¬ 
cies in this range. The most likely candidate is the n = 1 
toroidal mode; the detected frequen cy agrees ex tremely 
well with the most recent models dPiro 20051 . This 
would be fascinating in terms of energetics; the energy 
required to excite n > 0 modes is orders of magnitude 
larger than that required to excite an n = 0 mode, a 
testament to the strength of this particular flare. Other 
candidates include the cru stal spheroidal modes or the 
crust/core interface modes llMcDermott et al~llll988lh al¬ 
though the effect of a strong magnetic field on these 
modes remains to be calculated. If this mode is indeed 
a higher radial overtone, it could allow us to determine 
the depth at which the fracture occurred in the crust. 

The nature of the 18 and 26 Hz QPOs is less clear. 
The frequencies are too low for the fundamental crustal 
toroidal mode unless the neutron star parameters are 
very extreme. A torsional mo de of the core, restored b y 
the poloidal mag netic field ijThorri pson Dimca,nll2001h . 
was suggested bv llsrael et all (l2005jl as a candidate mech¬ 
anism for the 18 Hz QPO. Whether such modes could 
explain both of the QPOs is a matter for further study. 

The fact that the various QPOs are strongest at differ¬ 
ent times, suggests that we may be seeing two separate 
fracture/reconnection events. The first, associated with 
the main flare, excites the 626.5 Hz QPO. The second, as¬ 
sociated with the late time boost in unpulsed emission, 
excites the 92.5 Hz QPO. The idea that we might be 
seeing separate fracture sites is given additional weight 
by the observation that the two QPOs are strongest at 
different rotational phases. Such sequential rupturing is 
often observed in terrestrial seismology. It is interesting 
that the 18 and 26 Hz QPOs are strongest at the rota¬ 
tional phase where the 92.5 Hz QPO subsequently ap¬ 
pears. One might speculate that the magnetospheric os¬ 
cillations associated with the low frequency QPOs (which 
are triggered by the main flare) slowly weaken the crust 
at this point and trigger the second fracture. 

Understanding the means by which the modes cou¬ 
ple to the magnetic field, and hence modulate the X-ray 
lightcurve, is now critical. The fact that the modes ap¬ 
pear in different energy bands is intriguing, as is the fact 
that the 626.5 Hz QPO appears at a rotational phase 
when a large part of the surface should be obscured by 
the trapped fireball. Much theoretical effort is clearly 
still required, but prospects for neutron star asteroseis- 
mology are bright. 
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